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We use the Bodwin-Braaten-Lepage factorization scheme to separate the long- and short-distance factors that 
contribute to the decay rates of T, rib (S-wave) and Xb,hb (P-wave). The long distance matrix elements are calcu- 
lated on the lattice in the quenched approximation using a non-relativistic formulation of the b quark dynamics. 



In heavy quarkonium decays that involve 
quark-antiquark {QQ) annihilation, this annihi- 
lation occurs at short distances l/AIg). Bod- 
win, Braaten and Lepage ||l| have shown that this 
enables one to factor such decay rates into a sum 
of products of a short-distance parton-level decay 
rate with a long-distance matrix element between 
quarkonium states. The short distance pieces are 
calculated perturbatively, while the long distance 
parts are accessible to lattice calculations. To 
lowest non-trivial order in v"^, the square of the 
quark velocity, (w^ ^ .1 for bottomonium) 



where the X's represent states of light partons. 
The P's are the short-distance {p ^ Mq) parton- 
level decay rates. Gi, Hi and Hs are the long- 
distance {p ~ Mqv, E ~ Mqv^) matrix elements 
that we calculate on the lattice. 

In our lattice calculations we have used 149 in- 
dependent equilibrated quenched gauge configu- 
rations on a 16'^ X 32 lattice with (3 = 6.0. Heavy- 
quark, and hence quarkonium, propagators were 
calculated using the non-relativistic formulation 
of Lepage and collaborators |^ . We used the lat- 
tice version of the quark action that is based on 
the euclidean lagrangian 
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which is valid to the lowest non-trivial order in 
. We calculate the quark Green's function that 
obeys the evolution equation ||^ 

G(x, t + l)={l- i/o/2n)"[/t^t(l - ffo/2n)" 

xG(x,t) + Jx,o^t+i,o , (3) 

with G(x,i) = for t < 0, and Hq = 
-A(2)/2Mo - ho . Here A^^) is the gauge- 
covariant discrete laplacian, and Mq the bare 
quark mass, ho = 3(1 — uq)/Mo, where uq = 
(0|iTrC/p,a,|0)i 

The matrix elements we calculate are defined 
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CPHH^/2) D D i^\'P)/M^ (5) 

CP\^P^T-XX^T-i'\'P)/Ml , (6) 

On the lattice, we calculate the related quan- 
tities GI, HI, defined graphically below 





where the larger dots represent the "sources" , the 
small dot in the numerator is the appropriate 4- 
fermi operator, and the small dots in the denom- 
inator represent point "sinks" . For our calcula- 
tions we generate the retarded (Eqn. (3)) and ad- 
vanced quark propagators from noisy point and 
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noisy extended sources on each of the 32 time- 
shces. (This differs from our preUminary calcula- 
tions, in which the 4-fermi operator was used as 
a source.) Then as T, T' — > oo 



Gt(r,T') ^ Gi 
Hl{T,T') Hi 
Hl(T,T') 
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where Ris is the radial wave function of the IS 
state and R[p is the derivative of the radial wave 
function of the IP state. 

For bottomonium, we use input parameters de- 
termined by the NRQCD collaboration ||^ , which 
in our convention are: the bare b-quark mass, 
Mob = 1.5, the inverse lattice spacing, = 
2.4GeV, and the physical b-quark mass, Mb = 
2.06. In Fig. 0we show G\ as a function of T,T'. 
It is clearly very close to the vacuum saturation 
value of 1. In fact G| - 1 « 1.3 x IQ-^. HI dis- 
plays similar behaviour, but is more noisy. H^ 
is plotted in Fig. ||. We notice that it displays a 
fairly obvious plateau at small T, T', which de- 
generates into noise for larger values of T, T' . No 
improvement in H^ is obtained by using the ex- 
tended source. Fitting the plateau, we obtain 

Hs/Hi « 0.06 . 

This is somewhat smaller than the value obtained 
from a simple perturbative estimate However, 
this estimate comes from assuming that Hg, be- 
comes negligibly small when the momentum cut- 
off is Aqcd- If one assumes, instead, that Hs 
becomes negligible at a cutoff closer to the bot- 
tomonium binding energy, then the perturbative 
estimate is closer to the lattice measurement. Of 
course, the lattice-regulated Gi , Hi and Hg, differ 
from their continuum counterparts at 0{as), but 
since our methods are equivalent to using mean- 
field improved actions, these renormalizations are 
expected to be small. 

We have also considered the S-wave decays 
through next-to-leading order in v'^. To this or- 
der, Gi is no longer the same for T and r/b. How- 
ever, we would need an improved action in order 



to calculate these corrections. In addition, there 
is a second term in Eqn. (1), 

Fi{S)t'i{QQC'+^Sj^X) , (10) 

where T'l is another perturbative parton-level de- 
cay rate and Fi can be calculated on the lattice 
using the Lagrangian of Eqn. (2). In the vacuum 
saturation approximation, 

Fi{S) = (0|^t^|o)(o|V-t(^ D)\|0)/M^ . (11) 

On the lattice we measure Fj* , defined as 
M'Fi 



We find that 
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where non-covariant and covariant refer to 
whether we use ordinary derivatives (in coulomb 
gauge) or gauge-covariant derivatives in Eqn. 
(11). As with Gi, Hi and Hg, Fj* requires renor- 
malization. Fi mixes with Gi. Since Fi/Gi ^ u^, 
this mixing can be significant. We have calculated 
these mixings to 1-loop order. Preliminary esti- 
mates of the Fi's which take these mixings into 
account are 

F* {renormalized) — 0.76 — non-covariant (15) 
F* {r enor mail zed) — 0.62 — covariant. (16) 

Finally, in table we present some mass and 
wavefunction calculations which were incidental 
to our calculations of matrix elements. Clearly 
our numbers are inferior to those obtained by 
the NRQCD collaboration since we work only 
to lowest non-trivial order in v"^. However, they 
serve as a consistency check of our calculations. 

We are now in the process of repeating these 
calculations for the charmonium system at /3 = 
5.7 (/3 = 6.0 has too small a lattice spacing for 
NRQCD at the charmed-quark mass. Our earlier 
attempts [|| used charmed-quark masses which 
were too large.) The charmonium system affords 
the opportunity to confront our calculations with 
experiment, since there is already sufficient ex- 
perimental data to allow extraction of Hg,. In the 
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Table 1 

Properties of S- and P-wave bottomonium from our simulations. The lattice quantities include mean field 
renormalizations. The mass of the IS state is obtained by using M = 2{ZmMi, — Eq) + £'„ with Zm and 
Eq set at their mean field values. 





LATTICE 


EXPERIMENT 


Mis 

Mip - Mis 

\Ris{Q)? 

l^'ip(0)P 


9.2766(9) GeV 
0.434(9) GeV 
4.33(2) GeV^ 
0.75(7) GeV^ 


Mt = 9.46037(21) GeV 
M^, - Mr ^ 0.4398(7) GeV 
7.2(2) GeV^ 



future, we hope to extend these calculations to 
next order in v'^ and a^, and then to include the 
effects of light dynamical quarks. We are also 
calculating the complete renormalization matrix 
through 0{as) for the four operators discussed in 
this paper. 





Figure 1. as a function of T and T'. 
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Figure 2. iJ| as a function of T and T'. 
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